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Abstract 

We study the properties of the double-frequency sine-Gordon model in 
the vicinity of the Ising quantum phase transition displayed by this model. 
Using a mapping onto a generalised lattice quantum Ashkin- Teller model, 
we obtain critical and nearly-off-critical correlation functions of various 
operators. We discuss applications of the double-sine-Gordon model to 
one-dimensional physical systems, like spin chains in a staggered external 
field and interacting electrons in a staggered potential. 

PACS: 71.10.Pm; 71.10.Fd; 75.10.Jm 

Keywords: Fermions in reduced dimensions; Lattice fermion models; Quantized 
spin models. 

1 Introduction 

The problem of determining the asymptotic behaviour of a conformal field theory 
(CFT) under the action of a relevant operator is well studied and understood. 
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also in view of its relation to the physics of many quantum one-dimensional (ID) 
and classical two-dimensional (2D) models. A less studied case concerns a CFT 
perturbed by two relevant operators. Namely, consider a CFT subjected to such 
relevant perturbation, gO/^^, with scaling dimension Ag < 2, that turns it into 
a fully massive quantum field theory (QFT). Then add another relevant pertur- 
bation, XOaxj which, if acting alone, would also make our QFT fully massive. 
Although this is not as common situation as the previous one, we will show that 
it displays interesting features which may be realized in physical systems. 

Without loss of generality, we assume that < A^, < 2, i.e. that the second 
perturbation is more relevant then the first one. Moreover, we shall require, for 
the time being, that the operator product expansion (OPE) of these two operators 
is closed in the sense that it does not produce other relevant operators. The most 
general case will be discussed in a separate section. Then, the naive expectation, 
which prevailed until recently, is that no qualitative changes occur in this case 
with respect to the standard situation of a CFT perturbed by a single relevant 
operator: the low-energy behaviour of the theory would be governed by the most 
relevant operator and, at any rate, the theory would remain fully massive. This 
expectation is however not a general rule if the two operators exclude each other, 
that is, if the field configurations which minimise one perturbation term do not 
minimise the other. In this case, the interplay between the two competing relevant 
operators can produce a novel quantum phase transition between two massive 
QFT's through a critical (massless) point. It is intuitively clear that such a 
scenario is only conceivable if the ratio of the bare coupling constants, \g/X\, is 
large enough. 

Recently Delfino and Mussardo (DM) |l|] considered the double-frequency sine- 
Gordon (DSG) model, which is a Gaussian model of a scalar field $, perturbed by 
two relevant vertex operators with the ratio of their scaling dimensions Ag/Ax = 
4. The Hamiltonian density of the DSG model reads: 



4>R,L being the chiral components of the Bose field, $ = (j)R + (f)L- DM have shown 
that there exists a quantum critical line A = Xdg) where the DSG model displays 
an Ising criticality with central charge c = 1/2. 



-gcos(3<l>{x) - Asin[(/3/2)<l>(x)]. 



(1) 



where 



(2) 
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The purpose of this paper is to investigate in detail the critical properties 
of this transition and to discuss physical applications to spin chains and other 
quantum one-dimensional (ID) systems. While we shall concentrate on the DSG 
model in what follows, it should be noted that the phenomenon is more general - 
similar quantum phase transitions do happen in more complicated models; those 
shall be addressed elsewhere 

Apart from the practical interest in physical realizations of the DSG model, 
this problem is of relevance also from a pure theoretical point of view. Indeed, we 
have already mentioned that the Ising critical point separates two strong-coupling, 
massive phases; hence, by definition, its analytical description is outside the range 
of applicability of perturbation theory. In this paper we propose a nonperturba- 
tive scheme, which is essentially based on the mapping of the DSG model onto 
another equivalent model - a generalised quantum Ashkin- Teller model on a ID 
lattice, where the Ising critical point becomes accessible. 

The paper is organised as follows. In section 2 we briefly discuss the DSG 
model at the semi-classical level and outline our approach to tackle the Ising 
criticality in the vicinity of the decoupling point jS"^ = An. In sections 3-5 we 
present a complete description of the critical properties of the DSG model at 
the Ising transition. In section 3 we introduce a quantum lattice version of a 
generalised AT model which, apart from the conventional marginal inter-chain 
coupling, also includes a 'magnetic filed' type of interaction, hai(j2- We show 
that, in the continuum limit, this model is equivalent to the DSG model at 

~ 47r. In section 4 we employ a new (a-r) representation of the deformed AT 
model which enables us to correctly identify those degrees of freedom that become 
critical. In section 5 we use the results of the quantum-spin-chain mappings to 
determine the physical properties of the DSG model close to the criticality. 

In the next two sections we discuss applications of the DSG model to physical 
systems. In particular, we discuss Ising transitions in a spontaneously dimerized 
S=l/2 antiferromagnetic chain in a staggered magnetic field (section 6), and in 
the ID Hubbard model with a staggered potential (section 7). In section 8 we 
show that at (3"^ < 2n the second-order Ising transition may transform to a first- 
order one. Section 9 contains discussion of the results and conclusions. 

The paper is supplied with four Appendices which provide some details con- 
cerning bosonization of products of order and disorder operators for a system of 
two identical Ising models, the relationship between the quantum AT model and 
lattice fermions, a mean-field treatment of the a-r model, and calculation of the 
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correlation functions. 

2 The model and its quasi-classical analysis 

In this paper we will be dealing with the particular version of the DSG model 
corresponding to the case g > 0. We will also assume that A > 0, although its 
sign is unimportant due to the obvious symmetry $ ^ — $. 

For 2tt < f]"^ < Svr the Gaussian model is perturbed by two relevant vertex 
operators closed under the OPE |T|]. The existence of the Ising phase transition 
can be qualitatively understood via a quasi-classical analysis: inspecting the pro- 
file of the potential U{^) as a function of the ratio x = \/^g- At x = one has a 
periodic potential of a single-frequency sine-Gordon (SG) model. For x 7^ the 
period of the potential is doubled in such a way that in the region < x < 1 
it can be viewed as a sequence of double-well potentials with the local structure 
A^^ + 5$^, where A < 0, -B > 0. Precisely at x = 1 each local double-well po- 
tential transforms to a $^ one (A = 0), and (in the Ginzburg-Landau sense) this 
is a signature of the Ising criticality with the central charge c = 1/2. It should 
be stressed that the double-well structure of the potential W($) only occurs for 
(7 > 0. Hence the Ising critical point only exist for positive g. The case of (7 < is 
qualitatively different: here the A-perturbation removes the degeneracy between 
neighbouring minima of cos /5$ and thus leads to soliton confinement (similarly 
to the analysis in Ref|Q). The spectrum in this case always remains massive. 

Thus, for model (|I]) with (7 > 0, a plausible scenario is that the relevant pertur- 
bations naturally act on the two constituent Ising models of the starting Gaussian 
model, leaving one of them massless on the critical line A = \c{g)- (Quasi- 
classically \c{g) = 4:g- A better estimation which takes into account quantum 
fluctuations yields a power law Xc{g) ~ g'^, where u = (327r — / (IGvr — 
Indeed, DM have argued that the Ising transition is a universal property of the 
DSG model (0) as long as jS"^ < Svr. This makes it possible to consider the vicinity 
of the point = Att where the description of the transition greatly simplifies. 
Namely, it is well known 0, H that [3"^ = Att is the decoupling (or Luther- Emery) 
point of the ordinary SG model (A = 0), at which the latter is equivalent to a 
theory of free massive fermions. The scaling dimension 1/4 of the A-term in (|I]) 
indicates that this point is special for the DSG model as well, suggesting an Ising 
model interpretation. In this paper we shall be working in the vicinity of the 
decoupling point. Rescaling the field the DSG model can be written in an 



4 



equivalent form: 

Ti-DSG = V [{d^(t)Rf + {d:^(t>Lf \ - ld:^(t>R{x)d:^(t)L{x) 

- g'COS"\/47r$(x) - Asin v^$(a;), (3) 



where w = (fo/2)(/3/\/47r+\/47r//9) and 7 = {vq/ \/A7^l3){l3'^—An). It was originally 
noted by DM that this form of the DSG model can be related to a deformed 2D 
Ashkin-Teller (AT) model, which can be viewed as the standard AT model, i.e. a 
model of two marginally coupled 2D Ising models with order parameters ai and 
(72, extended to include a magnetic field type coupling, \a1a2. DM then argued 
that, in the strong-coupling limit (A 00), the system effectively reduces to 
a single Ising model which may become critical if the temperature is properly 
tuned to its critical value. This argument can be extended to finite values of A, 
provided that the two Ising copies are originally in a disordered phase, and thus 
explains the existence of a critical line in the t-\ plane {t = {T — Tc)/Tc > 0) 
along which the deformed AT model flows from the ultraviolet c — 1 fixed point 
to the infrared c— 1/2 Ising fixed point. 

We shall elaborate on the DM argument by explicitly constructing the map- 
pings between the DSG model and various spin models, with the aim to fully 
describe the Ising transition and, in particular, calculate correlation functions of 
the operators in the original DSG model. The novel points of our analysis are as 
follows. First, we concentrate on quantum lattice spin chain models rather than 
on their classical counterparts (in the transfer matrix sense), as this formulation 
profits from the use of the powerful apparatus of various spin operator transfor- 
mations like the duality transformation. Secondly, we identify the correct degrees 
of freedom in the deformed quantum AT model that become critical. Since the 
original Ising models enter symmetrically, this is a nontrivial step which is accom- 
plished via a 'change of basis' transformation. Let us denote the critical degrees 
of freedom by a and the remaining gaped degrees of freedom by r (see the main 
text for precise definitions). The strategy to calculate the correlation functions 
is then to express the DSG-operators in terms of the lattice a and r operators. 
At the Ising transition a and r operators asymptotically decouple, with the a- 
operators being critical and (some of) the r operators acquiring finite average 
values. This allows us to trace the relation between the original DSG-operators 
and those from the operator content of the underlying critical Ising model, which 
ultimately accomplishes a complete description of the critical properties at the 
Ising transition. 
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3 Relation between DSG model and deformed 
quantum Ashkin- Teller model 

3.1 Quantum Ising spin chain 

We start with recollecting some basic facts about the quantum Ising (QI) spin 
chain. The Hamiltonian of the QI chain describes a ID Ising model in a transverse 
magnetic field 0: 

^Q/M = -E(^«+i + A<)> (4) 

n 

where a" are the Pauli matrices associated with the lattice sites {n}. The Hamil- 
tonian Hqi defines the transfer matrix of the classical 2D Ising model 0. 

An important tool in studying ID spin lattice models is the Kramers- Wannier 
duality transformation that we shall make extensive use of in the sequel. Consider 
a dual lattice consisting of sites {n + 1/2}, defined as the centres of the links 
< ra, n + 1 > of the original lattice, and assign spin operators yUn+1/2 to the dual 
lattice sites. The duality transformation then relates the dual spins to the original 
ones as follows 

n 

the inverse transformation being 

n— 1 

'^n = n ^^j+l/2y = l^n-l/2f'n+l/2- {^) 

j=0 

In Hamiltonian (^, the parameters J and A are interchanged under the duality 
transformation, so that J = A is the self-duality point where the model displays 
an Ising criticality. 

The operators and /im+i/2' conventionally referred to as order and disor- 
der operators, are mutually nonlocal and play the role of 'string operators' with 
respect to each other. In particular, they commute for m < n but anticommute 
otherwise. These commutation properties make the duality construction a con- 
venient starting point for introducing lattice fermions. It is immediate to check 
that two objects 

Vn = ^n^n-\l2 — ^^n-l/2^ni i'^) 
Cn = iO"nAin+l/2 = ~^f^n+l/2'^n^ (8) 

satisfy the anticommutation relations for the real (Majorana) fermions on the 
lattice, with the normalisation r]'^ = = 1. Relations (0) and (||) are nothing but 



6 



the inverse of the Jordan- Wigner transformation, with the 'direct' transformation 
being of the form 

n-l 

< = ^CnVn, K = VnY[ [iCjVj) ■ (9) 

In terms of the Majorana fermions the QI Hamiltonian becomes 

^Q/ = iE [JUVn+l - Vn) - (A - JXnVn] • (10) 
n 

The next step is to take a continuum hmit. To this end one introduces a lattice 
spacing ao, treats x = na^ as a continuum variable, and replaces //„ and Cn by 
slowly varying Majorana fields, rj{x) and Ci.^)'- 

rjn V2a^r]{x), Cn yj^oCix). 

Notice that the factor \/2 ensures the correct continuum anticommutation rela- 
tions: {^7(3;), ?7(?/)} = {C{x)X{.y)} = ~ y)- The Hamiltonian density of ( PD[ ) 
then is 

Hqi = ivCd^r] - imC?7, 

with V = 2JaQ and m = 2(A — J). Performing a chiral rotation of the Majorana 
spinor 

. _ -^ + C . _v + C ...X 



or, inversely. 



transforms this Hamiltonian to a standard form: 

^ = '^i-^Rd.^R + ^lO^l) - im^R^L. (13) 

The relation between the QI model and the massive Majorana QFT is therefore 
summarised as follows 

hm Hqj= f dxH&)(x) (14) 

ao-+0 J 

3.2 Bosonization of the deformed quantum Ashkin— Teller 
model 

The standard quantum Ashkin- Teller (QAT) model is defined as a model of two 
identical QI spin chains, described by Hamiltonians HQi[ai\ and ifg/[cr2] [see 
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Eq.@], which are coupled via a self-dual inter-chain interaction: 

Hqat = HQi[ai]+ HQi[a2]+ H'^TWi,a2], (15) 

H'ATWl,a2\ = ^E(<n<n+l<n<n+l+<n^2,n)- (16) 

n 

We will be interested in a deformed version of this model which, apart from H'j^j,^ 
includes a 'magnetic field' type of coupling between the chains: 

Hdqat = Hqat - (^In'^ln (17) 

n 

We wish to establish a relationship between the deformed quantum Ashkin- 
Teller (DQAT) model (0) and the DSG model. These two models can be mapped 
onto each other using the Zuber-Itsykson trick||]. The idea is to associate the 
two QI chains with two copies of Majorana fermions, 

combine and into a single Dirac field and then bosonize the latter using the 
standard rules of Abelian bosonization. 

First we notice that, in the case of two QI spin chains, the Jordan- Wigner 
transformation (^-(§) should be slightly modified. This follows from the require- 
ment that the spin operators belonging to different chains should commute. While 
this is automatically true for the disorder operators fJ'2n+i/2 = 2;, y, 2;; a = 1, 2) 
because of their bosonic character, to ensure commutation between the order pa- 
rameters crf„ and cr|„ one has to introduce two anticommuting (Klein) factors 

{ki,K2} = 0, kI = kI = 1, (18) 

and replace (|^), (§) by 

Vi-,n = '*l'^l,nA''l,n-l/25 Cl,n = ^'^lC'"l,nA''l,n+l/25 

V2,n = '«2Cr2,n/^2,n-l/25 C2,n = ''''^2(^2,nf^2,n+l/2- (19) 

With the spin variables o"i^2 and /ii_2 subject to the duality relations (^, (|), 
definitions (|19D ensure the correct statistics for the Majorana fermions rja^n and 

Ca,n- 

In terms of the lattice Majorana fermions, the standard QAT Hamiltonian 
becomes 

Hqat = i 51 51 [JCa,n{Va,n+l - Va,n) - (A - J)Ca,nVa,n\ 
a=l,2 n 

+ ^ Cl,nC2,n ['7l,n+l'72,n+l + '7l,n^?2,n] , (20) 
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and admits a straightforward passage to the continuum hmit. The corresponding 
Hamiltonian density reduces to a theory of two interacting massive Majorana 
fermions 



nQAT{x) = J2 ^l?[ea(a;)] + 8Kao^^n{x)^2R{x)^iLix)U{x), (21) 



a=l,2 

which is obviously equivalent to the massive Thirring model for a single Dirac 
fermion. Using the standard rules of Abelian bosonization which are briefly sum- 
marised in Appendix A.l, one maps the QAT model onto a (3"^ = An quantum SG 
model with a marginal perturbation : 



Hqat ^ Ho [$] - 8Kaod^(f)Rd^(j)L - (m/ao) cos V47r$. (22) 

Notice that, the bosonization of the QAT model can be also achieved by first 
mapping the model onto a single chain of spinless fermions, and then bosonizing 
the latter. This alternative route is traced in Appendix |B[ 

Turning to the DQAT model, we observe that the h-term in (^) is nonlocal in 
terms of lattice Majorana operators. However, it is known that, for two identical 
Ising models close to criticality, products of two Ising operators belonging to 
different chains, (Ti(T2, ^1^2, o'i/i2, AiiO"2, all with scaling dimension 1/4, can 
be expressed locally in terms of bosonic vertex operators with the same scaling 
dimension (see, for instance, P] and references therein). In Appendix A. 2, we 
rederive this correspondence starting from the lattice theory of two quantum 
Ising chains and paying special attention to the Klein factors: 

Kn+l/2/^2,n+l/2 = : COS v^$(x) :, (23) 

^in(^2n = -iKiK2 : sm^/n^{x) (24) 
^inKn+i/2 = -i/ti : cosv/^e(x) :, (25) 
/^i,n+i/2^2n = i«^2 : slu v^e(x) : . (26) 

Here = —(pR + 0l is a scalar field dual to $. 

Thus, Eq. (^) establishes the continuum bosonized version of the "magnetic 
field" coupling term in ([T7|): 



hJ2'^in'^2n = i ( — ) ^I'^s / dx : sin\/7r$(a;) : . (27) 



Notice that algebra ([TSD of the Klein factors allows one to identify them with 
Pauli matrices, 

Ki = n, K2 = T2, K1K2 = m. 
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Therefore, in the continuum hmit, the DQAT model can be represented in a 
diagonal 2x2 matrix form: 



limHDQAT = dxn{x), 

ao->0 ^ ' 




with the two Hamiltonians 

- f— ) cosy47r$(x) =F [ — ) sinA/?$(x) , (29) 
Veto/ 

having the structure of the DSG model (^ and differing only in the sign of the 
coupling constant h. The identification of the parameters is as follows: 7 = 8KaQ, 
g = m/aQ, and X = h/ao (it is understood that when 0, K ^ 00, m,h 
so as to keep the DSG parameters finite). 

The 2x2 matrix form of Ti reflects the symmetry of the DQAT model with 
respect to the interchange of the two constituent quantum chains, (Vu), under 
which the DSG Hamiltonians 7i± transform to each other. Formally, this sym- 
metry appears as "unphysical" because the Hamiltonian Ti. commutes with T3. 
This allows one to set 

K1K2 = i, (30) 

and study the Ising transition by projecting the matrix DQAT model onto the 
subspace of a single DSG Hamiltonian 7i+ (the choice = — i, 'H —>■ Ti- 

would be as good as the above one). However, while is a conserved quantity, 
each Ki is not, implying that correlations functions of operators which contain the 
Klein factors may involve transitions between the two different sectors = ±1. 

We shall see below that some physical models are automatically mapped onto 
a single DSG model with a fixed sign of h (section 5), while some others reduce 
to the matrix form (|28|) (section 6). In the latter case, the effective Vu symmetry 
corresponds to a physical discrete Z2 symmetry of the model in question, which 
is spontaneously broken in the ground state, a typical example of this kind being 
a spontaneously dimerized phase. In such situations, the Hamiltonians H±, when 
considered separately, describe only excitations above each of the two degenerate 
vacua. On the other hand, the two-fold degeneracy of the ground state necessarily 
implies the existence of topological kinks which interpolate between neighbouring 
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degenerate minima of two different potentials l{±, corresponding to T-l±. However, 
in the entire parameter space of the model (PBD, including the Ising transition 
point, the sets of minima of the potentials U± always have a finite relative shift. 
This means that the topological Z2 kinks always remain massive and thus appear 
as irrelevant excitations as long as Ising criticality is concerned. 

Thus, when working with a single Hamiltonian 7i+, all operators representing 
strongly fluctuating fields at the Ising critical point are automatically taken into 
account. All "off-diagonal" operators proportional to ki or K2, such as those 
given by ( ]25|) and (p6D, can be dropped as they represent short-ranged fields at 



the critical point. 

4 Deformed quantum Ashkin- Teller model in the 
(cr, r) representation and the Ising transition 

The original two-chain representation of the Hamiltonian (|T7p is not the most 
appropriate one to describe the role of the coupling term P7|), but some general 
statements can already be drawn at this stage. First of all it is clear that, if 
the chains are originally in the ordered phase, the role of the /i-term would be 
exhausted by removing the degeneracy between the four ground states specified 
by the signs of the order parameters ai = and a2 = (o'fn)- The lowest-energy 
sector, determined by the condition haia2 > 0, will stay massive anyway. (For two 
ordered Ising copies, the h-teim gives rise to an effective longitudinal magnetic 
field applied to both Ising systems and thus keeping them massive.) Therefore, 
new effects can be only expected if the two chains are originally disordered. 

We can give an heuristic argument in favour of an Ising critical point that 
appears at a finite h, if at h = the chains are disordered. For the sake of clarity, 
let us make a duality transformation and replace the original model by a pair of 
two ordered Ising copies well below Tc, coupled by the interaction h^ifi2- Consider 
a single quantum Ising model (|^) at J A. In the leading approximation, its 
ground state is fully polarised. The disorder operator creates a domain 

wall by flipping all the spins within the interval 1 < j < n. At short distances, 
the domain walls behave like hard-core bosons which, in the dilute limit, are 
equivalent to spinless fermions with the dispersion 

e{k) = 2 J — 2A cos kaQ. 
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The coupling term /i/ii/i2 creates or destroys a pair of domains walls, one for each 
chain. Therefore, in the dilute limit, the effective Hamiltonian for the fermions 
describing domain walls is of the form: 

H = Y.Y. + h,ff J2 + {el k + C2,-fc) (31) 

k 1=1,2 k 

where hejj oc Ah/ J. The combinations of the creation and annihilation opera- 

=1. It is 



tors that appear in the coupling term reflect the fact that 
clear that the model (^) is a theory of four Majorana fermions, two of which 
are not affected by the coupling term. The remaining part of Hamiltonian can 
be straightforwardly diagonalized and reduced to two Majorana fields with the 
masses 

Ah 

m± = 2J ± const.—-. 

The Ising criticality is reached when one of the masses vanishes, which occurs 
when h ~ J^/A. In the dual representation, the above estimation of the critical 
field should be replaced by ~ A^/J. Strictly speaking, the domain walls do not 
posses a fermionic statistics but contain a Jordan- Wigner type phase, which can 
only be neglected in the extreme dilute limit. 

Although the above argument provides a qualitative explanation of the tran- 
sition, a satisfying description of the critical region can only be reached if we 
are able to properly identify the degrees of freedom which get critical and those 
which do not (and, of course, take into account the correct stastistics of the 
domain walls). For this purpose, let us introduce a new Ising variable 

< = <^ln^2,n- (32) 

Namely, let us switch from the original two-chain representation, with basic spin 
operators crf„ and cr|„, to a new one where the basic variables are and = (Tf„ 
(due to the Vu symmetry of the model, the choice and = cr|„ would be 
equivalent). In the original {cri-cr2) representation, the local (i.e. at a fixed lattice 
site n) Hilbert space of the two-chain model is spanned by the basis vectors |cri, cr2) 
which are eigenstates of and cr|: 

The new local basis \a, r) is defined as 

(T^|cr, r) = (T|cr, r), r^|cr, r) = r|cr, r). 
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where a = ai, t = aia2- Comparing matrix elements of the operators af^ and 
CT2„ in the two bases, we find the following correspondence: 



Z Z Z rr-^-r^ 



In. 

We also need to define the variables and Un dual to (T„ and r„, respectively. 
The pairs ((T„, fin) and (r„, z/„) should obey the duahty relations (||) and (^. Using 
these relations together with (|33D, one finds out how the dual spins transform 
under the change of basis: 

A''l,n+l/2 — f^n+l/2^n+l/2y f^2,n+l/2 — 

f^l,n+l/2 — f^n+l/2y f^2,n+l/2 — f^n+l/2^n+l/2- 

(34) 



In the a-T representation the DQAT model ([T7| ) transforms to another two- 
chain model which we call the a-r model: 

H[a,T]=H^ + Hr + H^r. (35) 

Here if^ is a QI Hamiltonian similar to (|^) but with different parameters: 

H. = Y: {-J«+i + Kal) . (36) 

n 

H-r is also of the QI type model but the magnetic field is nonzero both in the 
transverse and longitudinal directions: 

H^ = KY: - E + A<) . (37) 

n n 

Finally, the coupling term is of the Ashkin- Teller type: 

HaT = -Y. (j^n^^n+lTnTn+l + • (38) 



A mean-field approach to (|35|) is outlined in Appendix here we shall con- 
centrate on the large-/;, limit. For large h, the r-degrees of freedom freeze in a 
configuration where (r^) ~ 1, and (r^) ~ A//i. The cx-degrees of freedom are 
then described by an effective Ising model 



(39) 



which can indeed become critical when A ~ ^/j7i. Although this simple picture 
holds only when h ^ A ^ J, K , if a universal behaviour is to be expected, then 
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we are lead to conclude that, both at strong and weak h, the Ising transition 
essentially corresponds the situation when the a degrees of freedom go massless, 
while the r degrees of freedom remain frozen in a disordered configuration with 
both (r^) and (r^) nonzero. Notice that, with the above strong-coupling descrip- 
tion, we are unable to determine the critical line and even estimate the strength of 
the irrelevant operators close to this line. This means that, even though the uni- 
versal properties of the DSG model at the Ising transition, including the singular 
parts of physical quantities and critical exponents of the correlation functions, will 
be captured correctly, the nonuniversal parts, such as prefactors and subleading 
corrections to main asymptotics, are not to be trusted. 

Let us conclude this section by noticing two important facts which will prove 
useful in the next Section devoted to the correlation functions. 

(i) Quantum critical points are associated with gapless phases which are real- 
ized under certain conditions imposed on the parameters of the model. As 
soon as any of those parameters is shifted away from the criticality con- 
straint, the system becomes off-critical. The corresponding perturbations 
to the critical Hamiltonian are therefore relevant operators with confor- 
mal dimensions determined by the universality class of the critical model. 
Notice that there is no magnetic field coupled to in the exact lattice 
Hamiltonian (^5]). Therefore a departure from criticality by changing any 
one of the couplings {h, J, Aj, or K) will give rise to an Ising mass term. 

(ii) From Eg. , if we reasonably take \K\ < 2J, we arrive at the conclusion 
that, for h > h^, the a degrees of freedom are ordered, while they are 
disordered otherwise. 

5 Correlation functions 

In this Section we use the results of the above quantum-spin-chain type mappings 
to determine the physical properties of the DSG model close to the criticality. 

5.1 DSG operators in the (cr, r) representation: UV-IR 
transmutation 

In preceding sections we have shown that (i) the DQAT model can be mapped 
onto the DSG model which is a Gaussian free theory of the field in the 
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ultraviolet (UV) limit, and that (ii) the DQAT can also be mapped onto the 
{a, r) model which, at the Ising transition, essentially reduces to a single critical 
Ising model of the order field ct{x) and disorder field ^{x) in the infrared (IR) 
limit. Our aim here is to find out how the operators of the DSG model, originally 
defined in the vicinity of the UV fixed point, "transmute" when going from the 
UV limit to the IR limit. 



5.1.1 Current operators 

We start with holomorphic, or current, operators that are made up of additive 
(analytic and anti-analytic) chiral parts. Of physical interest are the vector and 
axial current densities which, in terms of the bosonic field of the DSG model, are 
defined as 

J{x) = jR{x) + JL{x) = ^dMx), (40) 



Ux) = jR{x)-jL{x) = -^dM^) = -^dM^)- (41) 

'TT ^/7^ 



In physical situations, J{x) determines the smooth part of the charge or spin 
density [J{x) = p{x)], while J^{x) describes the corresponding charge or spin 
current [J^i^x) = j{x)]. As follows from ( |1 14| ) , ( |115D , these can be expressed in 
terms of the Major ana fields and 

Jix) = i[^iRix)^2Rix)+^iLix)Uix)], (42) 
Mx) = i[U{x)U{x)-^iL{x)^2L{x)]. (43) 

Making the inverse chiral rotation from (^^,^£) to (?7",C")) {(^ = 1)2), we can 
define a local lattice operator 

-^n = ^ {'ninV2n + ClnC2n) , (44) 

which reproduces ( |1^ ) in the continuum limit. Using the inverse Jordan- Wigner 
relations (p!9| ) and transformations (|33D , (|3^) , we obtain: 

T Z Z ( Z Z Z 2^ 



1 



2^n (/^n+1/2 ~ f^n-1/2) (45) 

(here we have implemented our Klein factor convention (pOD). Using the fact 
that the r-field is noncritical and has a nonzero expectation value, we pass to 
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the continuum limit and thus find the expression for the current density at the 
infrared fixed point: 

J(x) (46) 

where C ~ (r^) is a nonuniversal number, and ^{x) is the Ising disorder field at 
the criticality. Thus, the UV-IR transmutation of the current density is given by 

f -^dMx) UV, , , 

[ Cd^fx{x) JR. 

Turning to the axial current J5, we notice that the latter is related to the 
vector current J via the continuity equation: 

dtJ{x,t) + dMx,t) = 0. (48) 

As a result, the IR form of J^ix) can be immediately recovered. Indeed, the 
Ising disorder field /i is a scalar field (with zero conformal spin) . This means that 
holomorphic properties of the vector and axial current operators are lost at the 
IR fixed point. Then the result (^61), together with the requirements of Lorentz 
invariance and continuity equation (^5[), leads to the correspondence 

f -^dMx) UV, , , 

\ -Cdtfiix) IR, 

with the same prefactor C as in Eq. (^Tj) , provided that the velocity is set f = 1. 

As an important consistency check, we still need to explicitly determine the 
lattice version of the current operator. One might naively conclude that, since 
under the inverse chiral rotation (P^ the axial current density transforms to 



J5(x) = -i [Ci(a;)?72(x) + r]i{x)C2{x)] , (50) 

it would be sufficient to replace the r.h.s. of ( ^0]) by its local (single-site) coun- 
terpart which one might identify with J^. This route is misleading because, in 
the lattice formulation, the particle current is always defined on a link (n, n + 1) 
and, therefore, should be determined from the equation of motion 

idtJn = [Jn, Hdqat] ■ (51) 

which, in the continuum limit, is supposed to reduce to Eq.(?8). It can easily 



be checked that the (density) operator, J„, commutes with all the interaction 
terms in the lattice DQAT Hamiltonian ([T7|) , since all of them, including the h- 
term, are made up of the density operators. Therefore the above commutator is 
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only contributed to by the Majorana kinetic energy term in (|17D and so is easily 
computed: 

^dtJn = Qn,n+1 ~ Qn-l,nj (52) 

where ^ 

Qn,n+1 = {Vl,n+lC2n + Cln'72,n+l) ■ (53) 

In terms of the lattice spins this reads 

In the IR continuum limit this becomes ~ dtfi since dtfi ~ i[fi^,HQi] ~ fi^, as 
immediately follows from the dual version of the QI spin chain Hamiltonian. Thus 
we arrive at the IR representation of the axial current given by (^9|) . 

5.1.2 Vertex operators 

The analysis of the vertex operators 

is in fact simpler than that of the holomorphic operators. There are two cases of 
physical interest: P = ±^/^r and P = ±v^47r. 
From (^, (|0D and (jS^) it follows that 

sin y/n^ ~ crio'2 ^ t I, (55) 

(/ being the identity operator), as the r model is always off-critical. This is a 
reasonable result since the above operator is directly present in the DSG Hamil- 
tonian (^. However, as we mentioned in the previous Section, the operator 
sin ^/^T^ corresponds to the departure from the criticality of the cr-model. So, 
this operator, though having a finite average value, should also possess an extra 
term which, at the critical point, represents a strongly fluctuating field (with a 
power-law decaying correlation function). A more correct version of the formula 
is therefore as follows: 

sin y/n^ ^ I + e, (56) 

where e is the energy density (or a Majorana mass bilinear) operator. This 
follows from the first observation (i) at the end of Section ^, which states that 
moving h from its critical value results in the stress-energy tensor renormalisation 
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Table 1: UV-IR transmutation of DSG-operators 



UV limit 


9.$ 






IR limit 






I + e 



and, more importantly, in the appearance of the Majorana mass, i.e. the Ising 
energy-density operator. 

Furthermore, (p3|) and the lattice fusion rule (|3^) give 

cos v^'^' ~ fJ'ifJ'2 ~ (57) 

Thus, the operator cosy^$ is the most divergent operator of the DSG model, 
with a nonzero expectation value at h < and vanishing upon approaching the 
Ising critical point as 

{fi)r^{K-hf\ (58) 

Finally, the behaviour of /? = iv^ivr operators is determined in full analogy 
to the above discussion: 

Our results on UV-IR DSG operators transmutations are summarised in table 

I. 

5.2 Correlation functions 

Having identified the operators at the IR fixed point, we can now determine 
leading asymptotics of the correlation functions. Most of the operators of physical 
interest can be expressed in terms of the Ising disorder operator /x at the Ising 
transition point and in its close vicinity. 

For physical applications of the DSG model, it is important to investigate the 
dynamical susceptibility defined as the frequency-momentum Fourier transform 
of the retarded auto-correlation function of the Ising disorder parameter: function 

roo /■OO 

D^^\uj,p) = -i dx dte-'P''+'^'{[n{x,t),fx{0,0)]). (60) 

J-oo Jo 

It is known that, at criticality, D(r) = Xjr^l^ where r = (r, x) (r = it, and v 
is set to 1). Furthermore, away from criticality, in the ordered phase, D(r) = 
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[Ai / n) Ko{mr) where Ai is related to the Glaisher constant. The crossover be- 
tween these two regimes, r <^ ,^ = 1/m and r > C,, respectively, is complicated and 
is described in terms of the Painleve theory If, following [^, we introduce, 

dlnD 

C = T 

dr ' 

then the function Q can be shown to satisfy 

which is related to Painleve V equation and can, in turn, be shown to produce the 
correct conformal and massive limits. As the exact expressions are uncomfortable 
for calculations, we shall estimate the Fourier transforms from the asymptotes of 
the correlation function D{r). Therefore, first we summarise known results on 
the limiting behaviour of this function. 
At criticality 

^o(r) = (61) 

It is understood that this asymptotics is still valid away from criticality in the 
region r <^^. 

In the ordered phase 

the large-r limit of which is 



D>(r) = — iro(mr), (62) 

TT 



D>{r) = ^^T^e— . (63) 



In the disordered phase, thanks to the work by McCoy, Wu, Tracy, and col- 



laborators |TD[, we know that 
Ai ( V 

D<{r) = — \m^r'^ K\{mr) — KQ{mr) — mrKo(mr)Ki{mr) 

+ lK^,{mr)\, (64) 



which is understood to be the connected part of the correlator (i.e. without the 
constant piece). The large-r asymptotics of this function can be obtained from 
the relevant expansions of the modified Bessel functions quoted e.g. in [0: 

D<{r) = ^,-,e-'-. (65) 

The retarded function (^) can be found by calculating the Fourier transforms 
of the above asymptotic forms with a subsequent analytic continuation to real 
frequencies. This procedure is outlined in Appendix 0. The results are as follows. 
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Criticality. Using ( |186|) and (|182D , one obtains 



n(R)( ^ _ 0I/4 r(7/8) [ 9{p^ - uj^) + cos(77r/8)g(^^ - p') 

^0 ^^.P) - ^ ^f(l78)1 |^2_p2|7/8 



2i sin(77r/8)e(cj2 

|l^2 _ p2|7/8 



(66) 



Ordered Phase. Here we have a particle pole (|181|) : 



A ' "2 A 

Df\uJ.p) = - ^ [5(c. + e,) + 5(u; - e,)] . (67) 



Disordered phase. Specialising to a = —1/2 in (|186|) we have: 



/(.;) = ^\uj^- el\e{ei - u;') + 2z^\uj^ - e^eiu;' - e^). (68) 

Therefore, as long as u"^ < p"^ + Arn^, the argument of the log-function in 
( |185D is real and positive, so the correlation function remains real without 
dissipation. Above the two-particle threshold (tu^ > + 4m^), the cor- 
relation function has a branch cut and the dissipative part will appear as 
follows: 

Df\co,p) = 

*ln| , U(p2 + 4m2-^2) 

' — uj^ + p'^ + Am?\ + 2m / 



Arctg ^""' P' ^!I}Le{u^ -p'-Am')]. (69) 
m J 

The excitations are therefore incoherent in this regime. 

6 Dimerized Heisenberg chain in a staggered 
magnetic field 

The DSG model exhibiting a nontrivial flow towards Ising criticality can be re- 
alized as an effective continuum theory for a number of quantum ID models 
of strongly correlated electrons, in particular quantum spin chains and ladders. 
In the context of spin systems, an effective DSG model can emerge within the 
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Abelian bosonization scheme when staggered fields breaking translational invari- 
ance, such as an exphcit dimerization (bond alternation) or a staggered mag- 
netic field, are added to an originally translationally invariant model with a gaped 
ground state[|. 

Perhaps the simplest example of this kind is given by the spin-1/2 Heisenberg 
chain with nearest-neighbour (Ji) and next-nearest-neighbour (J2) antiferromag- 
netic exchange interactions 

Hji~J2 = X! (^l^n ■ Sn+l + <^2S„ ■ S„+2) (70) 

n 

This model has been extensively studied during past years. If frustrating 
interaction J2 is small enough, the model maintains the critical properties of the 
unfrustrated Heisenberg chain (J2 = 0). At J2 > J2C — 0.24Ji, frustration gets 
relevant and drives the model to a massive phase characterised by spontaneously 
broken parity [p!3|, H. The ground state is dimerized and doubly degenerate. 



and there exist massive elementary excitations - topological Z2 kinks carrying 
the spin 1/2. At a special (Majumdar-Ghosh) point, J2 = 0.5Ji, the picture 
is particularly simple because the two Z2-degenerate ground states are given by 
matrix products of singlet dimers formed either on the lattice links < 2n, 2n + l > 
or < 2r2 — 1, 2r2 >. 

Under assumption that Ji ^ J2, the continuum limit of the Ji — J2 model 
can be considered, and the resulting quantum field theory is that of a critical 
SU{2)i Wess-Zumino-Novikov-Witten (WZNW) model perturbed a marginal 
current- current interaction 0]: 

271 

Hj,_j, = — {:3r-3r: + :3l-3l:)+i3r- Jl, (71) 

where 7 ~ J2C — J2 > 0. At 7 < 0, the perturbation is marginally irrelevant. 
However, at 7 > the effective interaction fiows to strong coupling, and the 

^While alternation of the nearest-neighbour exchange constants can originate from the spin- 
phonon coupHng, the case of a nonuniform magnetic field with a period 2ao (oq being the 
lattice constant) used to be regarded as unrealistic, not achievable in experimental conditions. 
Fortunately, the status of the staggered magnetic field has recently changed from exotic to 
legitimate. It has been shown that, in some quasi-lD antiferromagnetic compounds, such 
field can be realized as an intrinsic one. An effective staggered field can originate from the 
Neel ordering of one magnetic sublattice and being experienced by loosely connected magnetic 
chains which form another sublattice and remain disordered down to very low temperatures 

. A sign-alternating component of the magnetic field can be also effectively generated due 
to a staggered anisotropy of the gyromagnetic tensor, as it is the case for the Copper-Benzoate 
organic molecule [|12[. 
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system ends up in a spontaneously dimerized phase with a dynamically generated 
spectral gap[Q: mj_im ~ A exp(— 271/7), where A ~ Ji is the UV cutoff. 



Using Abelian bosonization, we can rewrite (ITT]) as a sine-Gordon (SG) model 



{d^^f + {d^Qf + ^d^^nd^^L - jTT^ cos V8n^. (72) 



The "hidden" SU(2) symmetry of this bosonic Hamiltonian is encoded in the 
robust structure of the last two terms in Eq.(^), parametrised by a single cou- 
pling constant 7. This fact enforces the SG model (|72D to occur either on the 
weak-coupling SU(2) separatrix of the Kosterliz-Thouless phase diagram (7 < 0), 
or on the strong-coupling SU(2) separatrix (7 > 0). In the latter case, quantum 
solitons with the mass nidim and topological charge Q = \ are identified with the 
Z2 dimerization kinks carrying the spin S = Q /2 = 1/2. In the remainder of this 
section we will be dealing with the massive, spontaneously dimerized phase. 
Consider the following deformation of the model: H = Hj^^j^ + H\ where 

H' = j2XaTr{Tag). (73) 

a=0 

Here g is the 2x2 WZNW matrix field with conformal dimensions (1/4, 1/4), 
and Ta are the Pauli matrices including the unit matrix tq = /. The scalar and 
vector parts of 



~ Tr (r^) ~ (^ cos v27r6, sin v27r6, — sin v27r$j , (74) 
eg ~ Tr (^) ~ cos -\/27r$, (75) 

constitute the staggered magnetisation and dimerization field of the S=l/2 Heisen- 
berg chain. Representing 

H' = Ae, + h, ■ n„ (76) 

let us consider the two terms in (^) separately. 

The role of weak explicit (spin-Peierls) dimerization in the spontaneously 
dimerized J1 — J2 spin- 1/2 chain has already been addressed by Affleck 0. The ef- 
fective double-frequency-sine-Gordon potential appearing in this case is different 
from the one studied in this paper (c.f. Eq.(|l])): 

Udimer = - J^~^ VStT^ + A COS VOjl^. (77) 

The A-term in (|7^) removes the degeneracy between the neighbouring minima of 
the unperturbed potential — cos y/Sn^ (i.e. between the two degenerate dimerized 
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ground states) and thus leads to confinement of tlie solitons. Tlie main pliysical 
effect is tlie spinon-magnon transmutation: deconfined spinous of the frustrated 
Heisenberg chain, carrying the spin S = 1/2, form bound states with S = and 
S = 1, the latter representing coherent triplet magnon excitations. 
Let us concentrate on the case of the staggered magnetic field, h^: 

H' = hs- n,. (78) 

Choosing = hsZ, we arrive at a bosonic model 

H = /Jo[$] + 7?^9A9,$L- 777^cosv^$ 

- /i,sinv^$, (79) 

in which we recognise the DSG model with the structure (|l]). From the analysis 
of the preceding section we conclude that the spontaneously dimerized chain in 
a staggered magnetic field has two phases separated by a quantum critical point 
at hs = h*. At hs < h* a "mixed" phase is realized, with coexisting dimerization 
(cs) 7^ and staggered magnetization (ris) 7^ 0. Notice that, as opposed to the 
case of uniform magnetic field that couples to the (conserved) total magnetisation, 
the dependence = ns{hs) shows no threshold in hs. Dimerization vanishes at 
the Ising critical point hs = h*s and remains zero in the "pure Neel" phase, 
hg > h*s. The critical field /i* can be estimated by comparing the dimerization 
gap mdim with the gap that would open up at 7 < 0: rrih ~ h'j,^^. So the critical 
staggered field is exponentially small: h^ ~ {rndimf^'^ ■ 

Since the spin SU(2) symmetry is broken by the (staggered) magnetic field, 
the total spin is not conserved, but the spin projection is. The latter cir- 
cumstance allows one to identify the spin of elementary excitations as the 
topological quantum number of the kinks interpolating between the nearest de- 
generate minima of the potential W($) in (^). According to the structure of 
W($), there will be "short" and "long" kinks, carrying the spin = | =F ^, 
where 6 = 6{hs) smoothly increases from 6 = at hs = to 6 = 1/2 at hg > h*. 
Therefore, in the mixed phase, the original massive S = 1/2 spinon splits into two 
topological excitations carrying fractional spins S^. These spins become 5*^ = 1 
and S*^ = at the Ising transition, and it is just the singlet kink which loses 
its topological charge and becomes massless at hs = /?-*. The existence of the 
fractional-spin excitations in the mixed phase [hg < h*) is nothing but the spin 
version of the charge fractionization of topological excitations found earlier in ID 
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Table 2: UV-IR transmutation of operators 



field 


UV 


IR 


uniform spin density 


J' ~ $ 




uniform spin current 


J ~ dt^ 


J ~ <9t/i 


dimerization 


Es ~ COS v^27r$ 


~ 


staggered magnetisation 


~ sin -\/27r$ 


~ / + e 



commensurate Peierls insulators with broken charge conjugation symmetry (e.g. 
cis-polyacetylene) ||15|| , and also in a recent study of a ID Mott insulator with 
alternating single-site energy [|16|. 

To estimate the behaviour of physical quantities at the transition, let us con- 
sider an anisotropic (7||,7_l) version of the model in which 

H = i7o[$] + ^9A5,$i--^cosv^$ 
27r [2-Kay 



- /i,sinV27r<l>, (80) 
The 7||-term in (|5D|) can be eliminated by an appropriate rescaling of the field, 

H Hq[(^] — — — — cos VSTTir* + hs sin V27rir<l>, (81) 

As already mentioned, universality arguments lead to the conclusion that the 
anisotropic model ( |8TD also incorporates the Ising criticality at some value of hg- 
Choosing -ft'(7|i) = 1/2, we reduce the perturbation to the form (^) 

H' = cos 74^$ + hs sin y/^^, (82) 

(27ra)^ 

discussed in detail in previous sections. Using the (cr-r) representation, we derive 
the relations in Table |^ describing the UV-IR transmutation of the physical fields. 
We have the following correspondence: 

hs < hi : disordered phase : (/i) 7^ 0; 
hs > hi : ordered phase : (/i) = 0. 

We see that dimerization is finite at hs < hi and vanishes as 

{es) (hi - hs)'/', 



24 



on approaching the critical point. The staggered magnetisation, on the other 
hand, is always finite in both phases. Its behaviour at the transition is determined 
by the subleading correction to the identity operator (see Table |I]): 

- (/^s-/^:)ln^^-^. (83) 

The logarithmic divergence of the staggered magnetic susceptibility at the tran- 
sition is similar to that of the specific heat of the Ising model: 

Next we consider the dynamical magnetic susceptibility. In analogy to (|60[), 
it is defined by 

/oo roo 
dx dte-'P''^'^'{[S'{x,t),S'iO,0)]) (85) 
-oo JO 

Using the continuum limit decomposition of the spin operators and the above 
glossary, we conclude that while the uniform magnetic susceptibility is readily 
given by 

x(cc^,g~0)~g2Z}(^)(^,?) (86) 

(the function D^^-^ has been extensively discussed in Section (|D), the staggered 
susceptibility is in turn related to the correlation function of the Ising energy- 
density operator. The latter object is a Majorana bilinear 

Therefore, calculating the staggered magnetic susceptibility reduces to a sim- 
ple task of computing the polarisation loop diagram for free, massive Majorana 
fermions. The result is 

, = ~ In (l - - 2^/ ^^arctg (^^-j— j j . (87) 

inside the gap (Icul < 2m) and 

X{u, vr) - In - l) - ^sign.. |vr + 2\/^^^^^arctg [jj^) | , (88) 

for \uj\ > 2m. 
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We notice that for h > hgit obviously is the staggered field operator = 27r) 
which alone dictates the physics of our effective DSG model (not only is it the most 
relevant operator, but it also has a large amplitude). It is therefore instructive 
to compare the results for the dynamical magnetic susceptibility of the DSG 
model we have found via Ising-type mappings with those for the SG model with 
a = 27r operator only, which were obtained in the paper Jl^ by means of the 



form-factor technique. As for h > hg we enter the ordered phase ((/i)=0), the 
uniform susceptibility (^,0) is of a coherent nature, very much in agreement 
with Jl^ ('magnon' contribution). The staggered susceptibility (|^,|^), on the 



other hand, is incoherent in excellent agreement with the kink-antikink continuum 
observed in Jl^ (the breather contribution is missed in our approach; it might 
be recovered as a bound state of Majorana fermions, but we shall not push our 
analysis beyond this point). No analogous comparison can be made for h < hg 
(no form-factor calculations are, to our knowledge, currently available for the 
= Svr SG model, and, even if they were, a comparison would have been of 
dubious validity, as the addition of the = 27r-operator qualitatively changes 
the spectrum from the start). 

7 Ising transition in the ID Hubbard model with 
alternating chemical potential 

The DSG model finds a number of interesting applications in the theory of ID 
strongly correlated electron systems. In this section we shall consider a partic- 
ular example of this kind which has been recently discussed in Ref. |]16[ - a ID 



repulsive Hubbard model at 1/2-filling with a sign-alternating single-site energy 
(i.e. staggered chemical potential). The Hamiltonian of this model reads: 

H = -tY. + h.c) + A E(-l)'ni,, (89) 



where Cjo- is the annihilation operator of an electron with the spin projection a, 
residing at the lattice site i, and rii^ = cl^Cia-. The model (|5^) was originally 



proposed in the context of quasi- ID organic materials |T8| and is also believed to 
be prototypical for ferroelectric perovskites [ll9| . 

In spite of its apparent simplicity, the model (|89|) reveals nontrivial physics. At 
f/ = it describes a band insulator (BI) with a spectral gap for all excitations. At 
A <^t, the low-energy spectrum of the BI is that of free massive Dirac fermions. 
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When the Hubbard interaction is switched on, the finite fermionic mass A makes 
the theory free of infrared divergences, so that the BI phase remains stable for 
small enough U. On the other hand, at A = 0, the Hamiltonian ( p9D coincides 
with the standard (translationally invariant) Hubbard model which is exactly 



solvable ||20| and is well known to describe a Mott insulator (MI) at any positive 
value of U, if the electron concentration n = {l/N)J2i,a^icr = 1 (the case of a 
1/2-filled energy band). The MI state has a finite mass gap rric in the charge 
sector induced by commensurability of the electron density with the underlying 
lattice. At energies well below m^, local charge fiuctuations are suppressed, and 
the low-energy spin dynamics of the model coincides with that of the spin-1/2 
Heisenberg antiferromagnetic chain, the latter possessing a gapless spectrum. At 
a finite U, the charge-gaped MI phase is stable against site alternation, provided 
that A is small enough |18| . 



Thus, the issue of interest is the nature of the crossover between the BI and MI 
regimes which is expected to occur in the strong-coupling region where the single- 
particle mass gap A becomes comparable with the MI charge gap rric. Starting 
out from the MI phase and decreasing f/ at a fixed A, one has to identify the 
mechanism for the mass generation in the spin sector. On the other hand, it is 
clear that the charge degrees of freedom should also be involved in the BI-MI 
crossover. Indeed, dividing the lattice into two sublattices, A and B, with the 
single- site electron energies A and —A, respectively, and considering electronic 
states of a diatomic (AB) unit cell in the limit t ^ [/, A, one finds a region U ~ 2A 
where two charge configurations, A^B^ and A^B^, become almost degenerate. 
This is the so-called mixed-valence regime where those excitations responsible 
for the charge redistribution among the two unit-cell configurations become soft. 
This means that, apart from the spin transition, a charge transition associated 
with vanishing of the charge gap at some value of U is also expected to occur. 



In Ref. ||16|| we have shown that the MI-to-BI crossover, taking place on de- 
creasing f/ at a fixed A, is realized as a sequence of two continuous transitions: 
a Berezinskii-Kosterlitz-Thouless (BKT) transition at U = Uc2 where a spin gap 
is djTiamically generated, and an Ising critical point at U = Ud < Uc2 where 
the charge gap vanishes. Assuming that f/, A ^ t, below we shall consider the 
effective low-energy field theory for the lattice model (^). We shall then briefiy 
comment on the spin transition and mostly concentrate on the Ising transition in 
the charge sector of the model which will be described in terms of a DSG model. 

The standard bosonization procedure (see e.g. p) allows one to represent the 
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Hamiltonian density as 

T^eff = Ti-c + Ti-s + "Hcs- 

Here the spin sector is described by the SU{2)i WZNW model with a marginally 
irrelevant current-current perturbation originating from the electron backscatter- 
ing processes ((7 ~ f/ao > 0): 

K = ^ (: Ji?^ ■ Ji? : + : Jl ■ Jl :) - 2g3n ■ Jl (90) 

where Jr^l are chiral components of the vector spin current satisfying the SUi (2) 
Kac-Moody algebra. This Hamiltonian accounts for the universal properties of 
the spin-1/2 antiferromagnetic Heisenberg chain in the scaling limit ||2l|. The 
charge degrees of freedom are represented by a sine-Gordon model for a scalar 
field 



rtc = — 
2 



mo 



Ui + {dM -^cos J87ri^,<l>, (91) 
J na * 



where uiq ~ g. The cosine perturbation is caused by the electron Umklapp pro- 
cesses (see e.g. 0). For a wide class of Hamiltonians with finite-range interactions 
(of which the Hubbard model is a member) the parameter < 1. This means 
that the model (|9l|) is in a strong-coupling regime, and the dynamically generated 
mass determines the gap rric in the charge sector. The charge and spin sectors 
are coupled by the A-term: 

2A 



Hcs = es sin j2nK,^, (92) 

na ^ 

where is the spin dimerization field of the S=l/2 Heisenberg chain, defined in 

Let us assume that we are in the MI regime with a gaped charge sector. Since 
the charge field $c is locked in one of degenerate minima of the periodic potential 
in (|9I|), 

(*c)^=J^m, me Zoo 



it follows immediately that, in the MI phase, the operator sin \^2n Kc^c appearing 
in (P^ ) does not represent a strongly fluctuating field; it is rather short-ranged 
at distances ~ Vc/ riic, and this explains the stability of the MI phase against a 
small A-perturbation. Assuming that the charge gap is nonzero, we integrate out 
the massive charge degrees of freedom to obtain the effective action in the spin 
sector. According to the unbroken SU(2) symmetry of the model, the effective 
Hamiltonian in the spin sector retains its form (0), but the current-current 
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coupling constant undergoes an additive renormalization. In the second order in 
A we obtain: 

9 9eff = g ~ C {/\/ m^f Vc 

where C ~ 1 is a nonuniversal numerical constant. The spin sector now resembles 
the Ji — J2 frustrated spin-1/2 chain (see section 5) with an effecive next-nearest- 
neighbour interaction J2 generated by the staggered chemical potential A. As 
long as Qeff > 0, the spin excitation spectrum remains gap less. However, decreas- 
ing f/ at a fixed A (or increasing A at a fixed U) eventually reverts the above 
inequality to g^jj < 0. In this case the current-current perturbation becomes 
marginally relevant, and a continuous (BKT) transition takes place to a sponta- 
neously dimerized insulating (SDI) phase with broken (site) parity and a finite 
mass gap in the spin sector. 

Thus, the condition g^jj = determines the spin transition point Uc2- Using 
the exact result for the charge gap in the small-f/ Hubbard modelpOd, 



we find that 



2nt 



1 + 



'lnln(t/A)\ 
. ln(t/A) j 



(93) 



In (t/U) 

The dimerization order parameter that becomes nonzero a.t U < Uc2 is defined 

as 

D = E(-l)* + h.c.) (94) 

and in the continuum limit its density is given by 



D{x) ~ cos ^J2nKc^cix)Trgs{x) (95) 

It is instructive to compare the parity properties of the SDI and BI phases. 
For models defined on a ID lattice, there are two parity transformations - the site 
parity (P5) and link parity (Pl) [^. The difference between P5 and Pl survives 
the continuum limit and shows up in two inequivalent parity transformations that 
keep the massless Dirac equation 

{dt + asd^) ipix) =0, ip = 

invariant. These are 

Ps : ^{x) ^ (Ti^{-x) (96) 
Pl : Hx) ^ (T2^{-x) (97) 
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Using bosonization rules for spin-1/2 Dirac fermions (see e.g. 0), one easily finds 
that both Ps and Pl lead to 



$c(x) ^ -$c(-a;), 

whereas the WZNW field Qs transforms differently: 

Ps ■■ 9s{x) ^ -gs{-x) (98) 
Pl : Ux) ^ 9s{-x) (99) 

Comparing ( P^D and (pSj), we see that the A-perturbation breaks but is in- 
variant under Ps (in fact, Ps is the symmetry of the Hamiltonian (^)), while the 
dimerization operator D breaks Ps but preserves Pl. This is easily understood 
by noticing that these two terms have the structure of two different fermionic 
mass bilinears, ip'^aiip and ilj'^a2ip, with opposite transformation properties with 
respect to Ps and Pl- Thus, the parity properties of the SDI and BI phases are 
different, and this is a strong indication that the passage from SDI to BI should 
be associated with a significant redistribution of the charge density (charge tran- 
sition) . 

In what follows, we will not be dealing with estimation of the transition point 
Uci ■ Referring the reader to Ref . [|TB| where the mechanism of the charge transition 



is discussed in the context of excitonic instability of the BI phase, here we simply 
claim that for the model (^), within the leading logarithmic accuracy, f/c2/f4i — 
1 = const/ In (t/ A), where the positive constant is of the order of unity. Let us 
instead focus on a nonperturbative description of the charge transition to the BI 
phase. Suppose we are in the SDI phase with both charge and spin sectors gaped. 
Adopting an Abelian bosonic representation for the SU {2)i WZNW model with 
a marginally relevant perturbation. 



2 



K + (5x$s) 



- -d^^.Rd^^sL + COS VS^'I's, (100) 

we can regard the effective Hamiltonian given by ( |100| ), (pT)) and (^2]) as a phe- 
nomenological Landau-Ginzburg energy functional, in the sense that all the cou- 
plings (mo, Kc, A, A) and velocities (f s, Vc) are effective ones obtained by integrat- 
ing out high-energy degrees of freedom. The effective potential is given by: 



W($c, *&s) = -/icCOsy87ri^c'^'c-/^sCOSv^$s-5siny27rfrc'^'cCOSA/27r$s (101) 
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with 

(J'c= — > 0, Hs = —, > 0, d = — . 



The robust ingredients of the potential ( |101|) are the vertex operators and the 
signs of the corresponding amphtudes. It can be shown that, as long as the 
parameter Kc is confined within the interval 1/2 < < 1, this potential, with all 
terms being strongly relevant perturbations, is indeed the most representative one 
for the model under discussion because no new relevant operators are generated 
in the course of renormalization. This is no longer true if Kc < 1/2, the regime 
which can be realized for an extended model with extra finite-range interactions, 
6-g- VJ2i''^i^i+i- 111 this case, new relevant vertex operators will be generated 
upon renormalization, and the continuous Ising transition in the charge sector 
can transform to a first-order one. We will comment on that in the end of this 
section. 

A simple analysis of the saddle points of the potential ^^(ec, £^s)|0 shows that 



the location of its minima in the spin sector. Eg = ^Tc/2n, and hence the spin 
quantum numbers of the topological excitations, are the same as in the BI phase 
{U < Uci). So the spin part of the spectrum in the SDI phase smoothly transforms 
to that of the BI phase. Therefore, being interested in the redistribution of the 
charge degrees of freedom in the vicinity of Ud, in (|101|) we can replace cos \/27r^s 



by its vacuum expectation value, (cos -\/27r$s) = ±co. The two signs here reflect 
the Z2 degeneracy of the dimerized ground state (spontaneously broken site par- 
ity; see Eq. (^)). Thus, the Hamiltonian of effective model describing the charge 
degrees of freedom can be represented in the following 2x2 matrix form: 

/ \ 

Hc;eff = : (102) 



where 



2 

- iXcCosJ9,TiKc^cThsm.J2'nKc^c (103) 



(here h = Scq). Notice that under P5 Hj^^^ ^ Hjr\ 

We have arrived at the matrix version of DSG model similar to 
To make contact with the DQAT model discussed in detail in sections 2 and 3, 
we shall consider Hc-gff in the vicinity of the point = 1/2. Setting = 
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1/2 (1 + 7o) and rescaling the fields 

transforms H^"^^ in ( |103|) to a form similar to Ti.± of Eq.(^). 

With this correspondence and all the results of the previous sections, we are 
now able to describe in detail the Ising transition in the charge sector. As already 
explained in section 2.2, this can be done by considering only the Hamiltonian 

When studying the physical properties of our model ( ^91) at the charge transi- 
tion point, it is important to remember that the disordered {h < he) and ordered 
{h > he) phases of the effective Ising model correspond to the SDI and BI phases 
of the electronic model (B^). First we consider the dimerization operator (pSj): 



P ~ cos A / 27: Kr^r COS V 27r$ , 



With the spin field locked in the SDI phase, (cos A/27r$s) = +cq and Kc = 1/2, 
this transforms to the charge polarisation field 

V ~ cos a/tt'^'c- 

In full agreement with the physical picture, from (0) it follows that the operator 
V, being the order parameter of the SDI phase, is the most strongly fiuctuating 
field at the Ising trasntion: 

P ~ /i (104) 

Its average value is nonzero in the SDI phase and vanishes as {V) ~ (hc — h)^^^ on 
approaching the charge transition point (remaining zero in the whole BI phase). 

Using (p5D, we find that the average value of the A-perturbation, which, at 
Kc = 1/2, is given by ~ sinvr^c is nonsingular across the transition and 
remains finite in both phases. 

The most interesting feature of the Ising transition is the UV-IR transmutation 
of the charge density Pc{x) and current Jc{x). At Kc = 1/2, the UV limit of our 
model represents a metallic state with central charge cuv = 2, described in terms 
of two massless Gaussian fields $c and $5. In this limit 

1 1 

pc{x) = -^d^^cix), Jc{x) = — ^d^Qcix). 



According to (^),(^), at the Ising criticality (c/^ = 1/2), 

Pc{x) = Cdxp{x), Jc{x) = -Cdtfi{x) (105) 
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With fi{x) representing the charge polarization field, Eqs.( |105| ) identify pc and 
Jc as the bound charge and polarisation-current densities, respectively. Such an 
identification is typical for insulators rather than metals. This could have been 
anticipated from the fact that a true metallic state with charge-carrying gapless 
excitations cannot be described by a single massless real (Majorana) field. 

An insulating (semi-metallic) behaviour of the model at the quantum critical 
point, reached at the charge transition, becomes manifest when one estimates the 
optical conductivity: 

a{uj) ~ -ujQmD^^\uj,0). (106) 

Here D^^\uj^q) is the retarded correlation function defined in (p(]|). Using the 
result (p6|), we find that, at zero temperature, aiu) displays a universal power-law 
behaviour: 

cTo(cu) ~ uj-^/^ (107) 

Although the optical conductivity is divergent in the zero-frequency limit, there 
is no Drude-peak contribution (~ typical of true metals. 

It is also possible to estimate the optical conductivity at finite temperatures 
(keeping in mind the a — r representation of the DQAT model, one should assume 
that T is much smaller than the mass gap in the decoupled r degrees of freedom). 
This can be done using conformal mapping from a cylinder [— oo<x<oo; 0< 
T < [3] onto a complex plane C . Starting with the T = asymptotics of the 
correlation function 

1 



{p{Zi,Zi)p{z2,Z2)) OC 



1/4' 



and using the above mentioned conformal mapping, one obtains 



(/i(x,r)/i(0,0)) OC 
It can be shown that I 



vrT 



sinh 7rT(x — it) sinh 7rT(x + it) 



1/8 



where 



5Jmx(c<j) 



p{x) 



2^7/4 



Ka-kt) 



IX 



(109) 



(110) 



^ (i6 

The final result for the temperature dependent optical conductivity at the Ising 
critical point is given by the following universal formula: 



<7{u, T) OC 



UJ 
2^7/4 



\AitT 



1111 
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At finite uj <^T, the frequency dependence of a{uj,T) is classical but with a 
quantum, temperature dependent prefactor: 



At u; ~ T, a{u!,T) reaches its maximum and then crosses over to its quantum- 



critical high-frequency {uj ^ T) asymptotics (|107|). 



8 First order transition at < 27r 

Until now we assumed that 27r < < Svr, so that no other relevant operators were 
generated upon renormalization. Indeed, already at < 327r/9, the operator 
sin [(3/3/2) which is generated by the OPE of the two operators present in 

(|I]), becomes relevant. However, such an operator does not modify the qualitative 
behaviour of the model, as one can easily realize by a quasi-classical analysis. On 
the contrary, the operator cos (2/5$(x)), which is also generated by the OPE, and 
which becomes relevant at < 27r, can modify the properties of the model in a 
relevant manner. Indeed, inspecting the quasi-classical potential 

U [$] = -g cos (3^ - X sin (^^^^ - V cos 2/5$, 

one finds that the Ising transition is turned by a sufficiently large V into a first 
order one. 

The capability of an apparently subleading operator to change a continuous 
transition to a first order one is indeed common to a variety of models. For 
instance, it is known that the critical line with non universal exponents separating 
the Charge Density Wave (CDW) and the Spin Density Wave (SDW) phases of the 
extended (U-V) Hubbard model, becomes a first order line at sufficiently strong 
coupling||2^. In fact, in the extended Hubbard model, the charge Luttinger liquid 
exponent Kc can be lower than 1/2, the point at which the second harmonics of 
the Umklapp scattering starts to be relevant. Therefore, at sufficiently strong 
interaction, it can indeed turn the CDW-SDW transition line into a first order 
one. A similar situation occurs with the charge transition in the electronic model, 
considered in section 7, when the latter is generalised to include a sufficiently 



strong nearest-neighbour repulsion 18 . 
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9 Conclusions 



In this paper, we have proposed a nonperturbative description of the Ising criti- 
cahty in the DSG model (|1|). Using the equivalence between the DSG model and 
a deformed quantum Ashkin- Teller model, valid in the vicinity of the decoupling 
point, 13"^ = An, we were able to identify the effective Ising degrees of freedom 
that asymptotically decouple from the rest of the spectrum and become critical 
in the infrared limit. This identification allowed us to describe the UV-IR "trans- 
mutation" of all physical fields of the DSG model and calculate the correlation 
functions at and close to the transition. We have also demonstrated the efficiency 
of our approach to describe Ising transitions in some physical realizations of the 
DSG model. 

We believe that our quantum-Ising-chain approach can be generalised to the 
case when the number of the constituent Ising models, coupled by the interaction 
hYljCTj, is larger than 2. Such situation can indeed be realized in certain SU(2)- 
invariant spin-ladders models which can be driven to criticality under the action 
of external staggered fields. In such systems, the quantum critical points may 
be not only of the Ising type but also correspond to SU(2)fc WZNW universality 
class. Such examples will be considered elsewhere 
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A Bosonization 



A.l Bosonization of Fermi fields, currents and mass bilin- 
ears 

We build up a Dirac field out of two Majorana fields, and and bosonize it: 



1- 


"6 + i6' 




[ V2 \ 



R,L 



-^exp(±iv^0fi,i) (112) 



To ensure anticommutation between the right and left components of the Fermi 
field, it is assumed that 

[<PrAl] = \. (113) 

Notice that a is an ultraviolet cutoff in the bosonic theory which actually appears 
as a short-distance regulator in the normal-mode expansion of bosonic fields. For 
this reason it does not need to coincide with the original lattice constant ag which 
appears in the continuum representation of fermionic fields. These two cutoffs 
are, however, related, as shown in Appendix A. 2. 

The chiral components of the U(l) current are defined as 

Jj, = ■.R^R:=i^,ji^^j,= ^d^<Pn (114) 



Jl = : LtL := i^iL^L = ^5,0L (115) 



Using ( [L 1 2| ) , one also finds that 

R^L = ■ 



l\/47r$ 



27ra 

^ (6r - (^iL + i^i) , (116) 



implying that 



cosV47r$ = ina {^iR^iL + ^2R^2l) (117) 
sin\/4^$ = -ina {^iR^2L + ^il^2r) (118) 

The DQAT Hamiltonian ([T7|) is symmetric under interchange (Vu) of the 
two chains. Let us find the transformation properties of all the fields under Vu- 
^From (|112|) it follows that interchanging the two chains leads to transformations 



R VR, VL ^ (119) 
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The chiral currents Jr^l^ and therefore the total {J = Jr + Jl) and axial ( J5 = 
Jr — Jl) currents change their signs: 



Jr(l) —Jr(l), J —J, J5 —J5 (120) 

Under V12 the scalar field ^ = (f)R + (pL and its dual counterpart 6 = —(f>R + (t>L 
transform as follows: 



e->-^-e (121) 



Notice that the symmetry properties of r.h. sides of ( [11 7]) and ( |118| ) are consistent 

with ([ml) . 

A. 2 Bosonization of products of order and disorder oper- 
ators 

Starting with the Jordan-Wigner transformation for two QI spin chains, which 
can be summarised as follows 

^an ~ KanVan, /^a,ra+l/2 ~ ~KanVa,n+l (122) 

(n \ n 

n ^aj Can, /^a,n+l/2 = H ^aj (123) 

f^L = i/^a < j Van, /^a,n+l/2 = ^aj^ ,Va,n+l(an (124) 

here we derive bosonized expressions for products of Ising fields belonging to 
different chains. 

• /^i,n+i/2 /^2,n+i/2 ■ Usiug the lattice definition ( |123| ), we have: 



n 



/^l,n+l/2/^2,n+l/2 " II ^lj'^2j (125) 



There are many ways to exponentiate the product in (T^), the naive assumption 
would be to use the identity 

and then replace in the exponential by iri(. This will bring us to a phase 
proportional to 

n 

iVlnCln + V2nC2n) 



37 



which, in the continuum hmit, reduces to 

Jo 



According to ( |117| ), the integrand represents a bosonic cosine operator, which 
is not what we would hke to have for practical purposes. To get a better rep- 
resentation, one has first to rearrange the four Majorana fields in the product 



a 



(126) 



It is readily seen that a product of two local Majorana fields, defined on the 
lattice, can be represented as 



?7i?72 = ±exp ( ±-r]ir]2 



TT 



(127) 



Therefore 



/^l,n+l/2/^2,n+l/2 " 11 " ^Xp 



128) 



Now we pass to the continuum limit by making chiral rotation ([T^) and using 



l^l,n+l/2f^2,n+l/2 



exp 



rx 

±7r / dy {^iR^2R + 6l6l) 
Jo 



rx 

exp =FiA/7r dy dy^{y) 



,=Flv^*(a;) 



(129) 



The l.h.s. of this relation is symmetric under 7^12, so must the r.h.side. As follows 
from (|121D , under Vu the field $ changes its sign. So the sign ambiguity in (|129|) 
is resolved by replacing the phase exponential by a cosine: 



/^tn+l/2/^2,n+l/2 = COSv^<l>( 



X) 



(130) 



<n+i/2 ^2,n+i/2- Usiug ([l23|) , wc havc: 



^ln'^2n 



{Klt^2) (ClnC2n) H ^Ij^ 



-Kll^2 ClnC2n (/^l,n+l/2/^2,n 



,n+l/2 



:i3ii 
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With the representation ( [1301 ) at hand, we only need to pass to the continuum 
limit in CinC2n and then make use of proper operator product expansions. We 
have: 



CinC: 



2n 



2aoCi(a;)C2(a;) 



(132) 



So 



lKlK20,o 



^ dx^{x) ^ sin V47r$(a;) 

IT Tca 



cos ^/^^(^{x + a) (133) 



Picking up the most relevant operators in the following OPE 

1 



dx^{x) : cos ^/^l^{x =F a) : 
: sin \/47r$(a;) :: cos A/7r'^'(a;) : 
one finds that 



: sin \pn^{x) : + 



sin ^^{x) : + ■ 



'ln^2n 



X] 



Intending to keep duality among (136) and (|130|) we set 



Oo 
TT 



(134) 
(135) 

(136) 
(137) 



and finally obtain: 

^In'^ln = -i'«i'«2 : sm^/^T^{x) : (138) 

Notice the important role of the Klein product kiK2- The l.h.s. of ( |138|) is 
Pi2-symmetric while in the r.h.side sin y^$(x) is antisymmetric (see ( |121| )). So 
the r.h.s. of ( |138|) is symmetric just due to the presence of K1K2. If we had 
replaced k,ik,2 by a constant, K1K2 — ^ i, we would obtain 

f^L^2n =: sin v^$(a;) : 

In this representation, the only way to ensure the Pi2-symmetry would be to 
impose a constraint that identifies $ and — $, in which case the scalar field $ 
would transform to an orbifold. 

• ^fn/^2,n+i/2- ^i^st wc represent af„/i^ ^+i/2 as 



f Z Z \ I z z \ 

,n+l/2 — \^lnl^l,n+l/2) \}''l,n+l/2l^2,n+l/2) ■ 
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The continuum limit for second product has aheady been found (see Eg. (|130|) ). 
The first product reduces to a Majorana field: o'infJ'i n+i/2 ~ "i'^iCin- In the 
continuum limit, with relation ( |137|) taken into account, 



Cln 



cos v^47r$/j(x) + cos \/A7r^L{x) 



(139) 



So 



COS v47r$i?,(x) + COS v47r$L(x) cosv^<l>(a; + a) (140) 



Making use of the following OPE 



: cos V 47r$/{(x) :: cos y/7r^{x + a) 



: cos v47r$L(x) :: cos ^/^l^{x + a) : 



- : cos ^/^TQ{x) : + 



we arrive at the result: 



^ln/^2,n+l/2 = "I'^l ■ COS v^e(x) : 



1411 



(142) 



/^i,n+i/2'^2n- Quite similarly one finds that 

f^ln+i/2^L = i«^2 : sinv^e(a;) 



(143) 



Using (|121|) , we see that under Vu the r.h. sides of Eqs. (|142|) and (|143|) indeed 
transform to each other. 
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B More on lattice fermions 



An alternative way to bosonize the QAT model, is to pass through a mapping 
onto spinless fermions, and bosonize the latter. 

We start by the quantum Ising model (]T0|). By means of the following unitary 
transformation 

Cn = -^{^Rn-^Ln), (144) 
Vn = -^iUn+^Ln), (145) 

the Hamiltonian (p!OD transforms onto 

= ^ ("2) ^ i^Rn+l — ^Rn-l) — ihn Hhn+l — ^Ln.-l)] 



+ X! (2^Ln — iln+l — C,Ln-l) — ^Ln {'^^Rn " ^Rn+l " ^fln-l)] 

n 

+ l{/^J-J)Y.iRniLn. (146) 

n 

We now consider the following quantum Ashkin-Teller Hamiltonian of two- 
coupled Ising chains 

Hqat = Hqi [fxi] + Hqi [a2] + H'^T [f^i , (^2] , (147) 

where -ffg/'s are Ising Hamiltonians [see Eq. ( |14(j| )] for two spin species, ai and 
(T2, and 

H'aT = KY: (<n<n+1^2,n^2,n+l + <n^2,n) , (148) 

n 

is the coupling term. For each spin species we introduce Majorana's fermions. In 
order to make ui and 02 commute we multiply each Majorana's fermion for the 
chain 1 by another Majorana ki, and for chain 2 by k,2- That is, if a = 1,2, we 
define 

Ca,n = l^a (i<„/i:,„+l) , ?7a,n = l^a (<„/i:,„_l) , (149) 

as well as their R, L components 

Cl,n = -^{^Rn-^Ln), (150) 

Vl,n = -^{Un + ^Ln) , (151) 

C2,n = iVRn - VLn) , (152) 

?72,n = iVRn + VLn) ■ (153) 
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Let us concentrate for the moment onto the two Ising Hamiltonians, which are 
bihnear in the Majorana fermions. We notice the following general property 
(p, q = R,L) 

Cpn^qm ~l~ ^pn^qm [('Cpn ^^pn) i^^qm ~l~ i^qm) ~\~ i^^pn ~\~ i^pn) i^^qra i^qrn)\ 



^pn^qm ' ^pn^qm 



where we define the Fermi operators 

_ 1 

CR{L)n - 2 

Therefore, through ( |146| ), we find 
Hqi [ai] + Hqi [a^l = 



^R{L)n + iVR{L)r 



(154) 



+^ r, X! \^^Rn {^(^Ln '^Ln-l) ^^Ln i^^^Rn %n+l %n~l) H .C. 



+2i (Aj -J)J2 CruClu - CluC 



Rn- 



(155) 



The analogy with a fermionic model on a lattice is already apparent. Let us 
make this analogy more firm. We consider the following spinless fermion model 
on a lattice of 2L sites 

2L-1 

H = -t 4cn+i + H.c. 

n=0 
L-1 

n=0 



2n+l 



We make the following unitary transformation 

(-1) 



C2n+1 



V2 

til 
V2 



{cRn + CLt. 



[CRn - CLt^ 



where the right- and left-moving fermions are defined on a lattice of half the 
number of sites, i.e. L. The Hamiltonian becomes 
t 



H 



^Rn+l ^Rn-1 



c\„ I C 



^Ln (2Cr„ ^Rn-l 



(156) 
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Let us add to this Hamihonian a staggered hopping 

6H = -Aj2cLc2n+i + H.c. (157) 

n 

= zAY^ci^c^^-H.c. (158) 

n 

We notice that (|156|) plus ([158D coincide with (|I5|) if 2 J = t and 2(Aj - J) = A. 

Now we focus on the coupling term ( |148[ ) , which, in terms of Majorana's 
fermions reads 

H'aT = -^Xl ClnC2n {Vln+lV2n+l + ^ln'72n) • 
n 

Transforming into right- and left-moving Majorana fermions we find that 

ClnC2n = - {^Rn — C,Ln) {VBn — VLu) 
= -i {PRn + PLn-l- A„) , 

where we used ( |154|) , and we defined Pp„ = c]^„Cp„ (p = R,L) as well as A„ = 
c^HnCin + H.c. We can equivalently show that 

VlnV2n = -i {pRn + PLn " 1 + A„) . 

Therefore the coupling term can be written as 

H'^T = H {pRn + pLn-^- A„) (pnn + pLn " 1 + A„) 

n 

- (pjj„ + pl„ - 1 - A„) (pr„+i +pi„,+i - 1 + A„+i) . (159) 

n 

On the other hand let us consider the simplest nearest-neighbour interaction for 
spinless fermions, namely 

^ '^2n'^2w'^2ra+ 1^2^+1 + (^2n+l(^2n+l(^2n+2(^2n+2- 
n 

This interaction term turns out to be equivalent to ( |159| ), apart from a chemical 
potential term, provided 4:K = —V . Notice that, if one considers an Ashkin- Teller 
coupling which is not self-dual, e.g. 

H'aT = KlY^ (^ln(^ln+l(^2,nO-2,n+l + ^2^2 ^l,n^2,n> 
n n 

this translates into a staggered interaction for the spinless fermion model. This 
interaction, even without an explicit dimerization, is able to gap the fermionic 
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spectrum, thus showing the importance of self-duahty, even at the level of the 
coupling term ( p.47| ), to get a critical behaviour. 

Therefore we have shown the equivalence between a model of two coupled 
Ising chains in a transverse field, given by the Hamiltonian ( |147| ), and a model of 
spinless fermion with nearest neighbour interaction and dimerized hopping. 

An interesting point in this respect regards quantum numbers. The spinless 
fermion model has for instance conserved number of particle. In p is the density, 
in the reduced chain we must have that 

\^PRn + pLn = '2^P, (160) 
n 

is conserved. On the other hand 

7 Y.{.PRn + pLn - 1) = 2p - 1 
^ n 

= ~777 X! ClnC2n + VlnV^n, 
^-^ n 

which, in terms of Ising variables, implies the conservation of a very non-local 
operator, which includes the spins and their dual counterparts. 

B.l Operator identities 

We can build up several operator identities in the lattice representation. 

(1) The simplest one is obtained by the identity 

2 {pRn + PLn - = i (ClnC2n + ^ln6n) , (161) 

which relates the density operator of the spinless fermions to a particular 
bilinear of Majorana's fermions. 

(2) The other identity derives from the equality 

2 {cRnCLn + H .c}j = 1 (6n6n " ClnC2n) , (162) 

which relates the charge density wave operator of the spinless fermions to 
the Majorana's fermions. 

(3) The last identity is obtained by the equality 

- 2i {cRuClu - H.C.^ = i (Cln6n + C2n6n) , (163) 

which provides a relation between the dimerization operator and the Majo- 
rana's fermions. 
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C Mean-field treatment of the (cr, r)-model 



In oder to devise an appropriate mean-field scheme, let us formally rewrite Hamil- 
tonian flSBf) as 



H[a, t] = H,nf [a] + Hmf [r] + i//i„c[fT, r] + const. (164) 
Here the two terms of the mean-field Hamiltonian are: 

Hmf[a] = - E (-/X^+i + AX) (165) 

n 

Hmflr] = - E + K + Arr:) (166) 

n 

where 

J. = J (l + (r^„%i)) , J. = J(«+i) - K 

A, = A{r:)-K, A, = A(l + «)) (167) 

The coupling between the fiuctuations, is given by 

n 

+ A«-«))(r:-(r:))] (168) 



In the leading approximation, fiuctuations described by (|168| ) are neglected 



and the a and r degrees of freedom decouple. We notice that the parameter h 
appears only in Hms\r\ where it plays the role of an effective longitudinal field. 
Therefore the model ( |166| ) is a quantum counterpart of a 2D Ising model in a 
nonzero external field which is known to be always noncritical. Thus, as expected, 
the r-chain is always gapped. The condition for the cr-chain to be critical then 
reads: 

+ =A(r:)-ir (169) 

Despite the fact that the Ising model in the field is supposed to be integrable 
in the continuum limit (see references in |1^), the explicit expectation values 
of the r operators appearing in (|169|) are not known. To make further analytic 
progress possible, we choose a special value of the field rescaling parameter K, 

at which = 0, and interaction between the neighbouring r-spins in ifmj[r] 
vanishes. This trivialises the r-model and leads to 

(0= , ^ , (r-) = , (171) 
+ A2 Jh^ + A2 
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Self-consistency requires the knowledge of the averages related to the cr-model. 
Specialising to the critical point and making use of the known results H, 



Z Z \ / X\ 



IT 



we obtain 



Then the condition ( |169D reduces to 



i + {Tr = {A/j){T^)-- 



TT 



(172) 
(173) 

(174) 



This equation determines the critical value of /i at a given ratio A/ J. Introducing 
the quantity 

VWTc^ ^1,2 

X = ; , (7=1 + - 

J TT 

we obtain a quadratic equation 

(C + l)x^ - C{A/Jfx + C2(A/J)2 = 0, 

whose solution determines the critical value of h: 
( hr 



1 



CA) 4(C + 1) 



-) +4(C 



(175) 



The requirement that the r.h.side of ( |175|) is positive yields a restriction upon A: 

(176) 

Recall that we already assumed that at /i = the two QI chains of the QAT model 



A 2 
->1 + - 

J TT 



are disordered: A > J. The restriction (|176|) tells us that, for the Ising criticality 
to be reached at some critical value of h, the original disordered ai and a2 chains 
should be far enough from their original critical point. In fact, as follows from 
( |176| ), the value of the Majorana mass, when estimated as m ~ 2(A — J) ~ 4 J/ tt, 
turns out to be of the order of the ultraviolet cutoff. This estimation reflects the 
already mentioned fact that the DM transition is indeed not a weak-coupling one 
from the standpoint of the DQAT model. 

One could, in principle, perturbatively calculate the fluctuation corrections, 
i.e. those originating from ( |168| ), to the criticality condition. As the model is 



only tractable under fine tuning (|170| ) of the parameter that is not expected 
to generalise the qualitative picture obtained in the mean-field approximation. 
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D Correlation functions 



Let us start calculating the Fourier transform of D^^\uj,p), given by (§3), with 
the remark that using the Bessel functions (despite some convenient Fourier trans- 
forms) is beyond the accuracy. The analytic properties (and the issue of coherent 
particle poles) should not depend on fine details of the correlators at r ~ ^ = 1/m. 
Therefore what one ought to do is to study analytic properties of the integral 

oo 

/-(g) = J r"-^/2^rJo(gr)e-"^'' (177) 



in terms of which the Fourier transforms, 

1 °° 
D{q) = -J d^fe''^D{r) = nj rdrJo{qr)D{r) , 



q= {ujn,p), q = ^JuJl+p'^, are given by 

In the area of convergence, we have 

= ^^^;;T^^("/2 + 1/4, «/2 + 3/4; 1; -q'/m') (179) 

where F{a, b; c; z) stands for the hypergeometric function. The character of the 
singularity at z = 1 can be determined by using appropriate transformation for- 



mulas for the hypergeometric functions and the F-function doubling formula [25 



T^(nM (2m)"-V^r(a) 1 

where only the leading, most divergent term is retained. 

This calculation explains how the correlation function has a branch cut at the 
threshold, unless a = 1 which, of course, corresponds to the disordered case: 



7T A 

D>{q) = (181) 



where the coefficient, as expected, is the same as one finds by Fourier transforming 

Knirar) = — / d^o- 



' q'^ + m'^ 
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Although this calculation is instructive, formula (|18CI|) doesn't solve all our 
problems. So, the critical correlation function has nothing to do with the m —>■ 
limit of this formula, as the main contribution to the integral comes from a 
different spatial domain. Therefore we should rather return to the general formula 
( |179| ), set m = there, and analytically continue beyond the convergence domain 
(-1/2 <a< 1/2). The result is 

The correlation function in the disordered phase is also outside the validity 
of ( 11801 ) and, furthermore, the integral ( |179|) is logarithmically divergent at the 
lower limit in this case. It should therefore be regularised. One way to do it is as 
follows H: 



■2mr 



uq" ^ V ? J \ J \ J 

This regular isat ion is, however, not good enough. That is because the integral 
still diverges in the q ^ limit as 



In 



V4m 



while the physical regularisation is finite in the q ^ limit and expands as: 

- />io(,r)e-- ^ / ^e-' ^f jg, L^)' (184) 



Since different limiting processes must lead to the same result, one can simply 
substruct the unphysical g — > divergence from ( |183| ). This leads to the following 
result for the correlation function (up to an unessential additive constant): 

, , 35v4i / 4m \ , , 

D<(q) = n ^==^ 185 

Turning to the problem of analytic continuation, we recall that, according to 
the standard approach 
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in the upper half-plane. Put another way, one has to substitute 



(that is provided the resulting function has no 'accidental' poles and is otherwise 
regular in the upper half-plane - there is no general solution to the analytic 
continuation problem.) The resulting function is, by construction, analytic in the 
upper half-plane and automatically satisfies the Kramers-Kronig relation. 

We need to analytically continue formula ( |180D , or equivalently the function 

1 



It is easy to see that 



ejel - uj^) + cos(7ra)g(^^ - 



,2 _ c2|o 



2ism(TTa)9(iu — e:,, 

^ ^ ^ (186) 
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